Low-energy t^N interactions: scattering lengths and resonance parameters 
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INTRODUCTION 



^ ' While most of our knowledge of the N and A baryons has come from ttN elastic scattering and 

, photoproduction, the TyN channel has been crucial in determinations of the iV(1535) properties. In 

ttN elastic scattering, this resonance signal is masked by a sharp cusp due to the opening ?7N channel, 
while in the reactions 7r~p — > rjn and 7p rjp, the A'^(1535) is associated with a rapidly increasing 
cross section near threshold. As a result, the incorporation of eta-production data in multi-channel 
I fits has allowed more reliable determinations of the N{1535) resonance parameters. 

f"^ . The nearby iV(1520) resonance, while not strongly coupled to the 77N channel, gives an important 

contribution to some eta-production observables through interference effects. The Particle Data 
Group P] estimates the ratio T^n/^ tot to be 0.0023±0.0004, a value determined mainly by the 
multi-channel analysis of Penner and Mosel 0|- In Tr~p rjn, the effect of this resonance is 
O \ visible in the departure from purely S-wave behavior with increasing energy. While the iV(1520) 

2 ' contribution is small, the effect is magnified through S-D wave interference with the dominant 

iV(1535) contribution. In rjN photoproduction, this interference effect is particularly evident in 
measurements with a polarized beam (S) 

The 77N interaction has also been studied extensively due to a strong attraction at low energies, 
^ observed in most analyses, which could possibly lead to the existence of bound state 7y-mesic nuclei 

' ^ IS ' Unfortunately, the rjN scattering length cannot be measured directly. Instead, cross sections 

for 72_incson production, 7r~p rjn and 7p rjp, have been studied [a S H IS Qi ^1 Qi 
3 16, 17, 18, 19, In Ref. j21j|, it was demonstrated that the real part of the scattering 
length cannot be extracted directly from the low-ener gy e ta- p roduction cross sections and a model 

analysis is needed. Previous analyses li IM la |2ll, MllJll 111 111 111 IM IM IIS IM 111 IM 



1421 l43l |4J, |45|, kg, l47l l48l l49l |50j have found a large spread for the real 
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part, from a negative value |2l| to 1 fm [50j (Table QJ. Reasons for the spread include the rather 
old and conflicting 77 meson production data and differing TriV elastic scattering amplitudes from the 
Karlsruhe [sj and SAID |53| groups. However, the largest factor appears to be the model used in 
analyzing the data. 

In recent years, the eta-production database has seen significant improvements. In this paper, 
we present the results from analyses of these new data. The most recent data for the reaction 
7T~p rjn are described in Section |n] In Section ITTll we described the formalism associated with 
our fits. Results and comparisons with previous determinations are tabulated in Section llVl Finally, 
m Section El we summarize our findings and consider the open questions which will require further 
work. 
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II. THE EXPERIMENTAL DATABASE 



Until recently, the Tr~p — » rjn cross section database contained mainly old and often con- 
flicting measurements (see Ref. for details), with no polarized measurements existing below 
1.1 GeV/c jH^]. These older data have been reviewed by Clajus and Nefkens [S^. More recently, 
differential and total cross sections for Tr~p rjn near threshold have been measured using the BNL- 
AGS Facility (E909) 0. Data were obtained from threshold {p^ = 684.5 MeV/c) to -770 MeV/c. 

The most recent ■K~p rjn experiment was also performed at BNL, but in this case, using 
the Crystal Ball spectrometer (now moved to MAMI at Mainz). Cross sections were measured from 
threshold to 747 MeV/c (E913/914) 0. The total and differential cross sections from these two BNL 
measurements are compared in Fig. ^ Normalization issues still remain for the differential cross 
sections. However, the Crystal Ball distributions show a much smoother variation and a clear onset of 
higher partial-wave contributions with increasing energy. This feature, which is a vast improvement 
over previous measurements, allows an improved separation of the iV(1520) contribution. 

Traditionally, the total cross-section is plotted as a function of the pion laboratory energy 
[Fig^h)]. This view shows the sharp growth above threshold which is usually attributed to the 
dominance of the iV(1535) resonance, having a mass close to the rj production threshold (a/s = 
1487 McV) and a strong coupling to the rjN system. In Fig. |21 we instead plot the total cross- 
section as a function of the rj cm momentum p*. From the figure, we see that data can be described 
very well by a linear fit (dashed line) . This is due to the 5- wave dominance of the total cross section. 
From the slope of the best-fit line, a restriction on the imaginary part of the rjN elastic scattering 
amplitude A^jq can be found. 



The optical theorem leads to 



As a result, we have 



Pr] 

3 P^^ P*] 



3p* ^ 

ImA,,N > -^-^ CT^-p^,,n. (2) 
Snp* 

Using a linear fit, the recent E909 threshold data give 

\^.-p^nn = 15-2 ± 0.8 jih/MeV 

p* 

ImAriN > 0.172 ± 0.009 fm, (3) 
which can be compared with a previous output from Ref. Q 

—<y-K-v^vn = 21.2 ± 1.8 /ib/MeV and 

ImAriN > 0.24 ± 0.02 fm. (4) 

It is commonly believed that the iV(1535) resonance dominates the rj production cross section. 
This resonance mechanism results in the imaginary part of the rj pion-production and the rjN elastic 
scattering amplitudes being determined mainly by the A'^(1535) resonance parameters. But this is not 
the case for the real part. The real part of the resonance amplitude goes to zero at the resonance 
position. Therefore, the real part of the rjN scattering length strongly depends on nonresonant 
processes. For this reason, a multichannel analysis is favored in determining the r]N scattering 
length. Our approach is described in the next section. 
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III. COMBINED ANALYSIS OF ttN ELASTIC AND tt p ^ r/n DATA 

Our energy-dependent partial-wave fits are parametrized in terms of a coupled-channel Chew- 
Mandelstam K matrix, as described in Ref . [s^ . This choice determines the way we modify energy 
dependence and account for unitarity in our fits. Data for TriV elastic scattering have been fitted 
up to 2.1 GeV in the pion lab kinetic energy. Data for the reaction Tr~p -qn have been included 
from threshold up to 0.8 GeV. Constraint data have also been included, in order to ensure that 
the resulting fit produces elastic TriV amplitudes satisfying a set of forward and fixed-t dispersion- 
relations. This fit to data plus constraints must be iterated until a stable result is obtained [53l |. 
Finally, we have included ttA and pN channels to account for unitarity, but have not explicitly fitted 
data of this type. 

The iV(1535) resonance couples mainly to ttN and rjN ^ with a much smaller branching fraction 
to ttttN, and our results were not sensitive to the choice of additional channels. For the iV(1520), 
however, there is a substantial inelastic branching to ttttN, split mainly between pN and ttA. We 
therefore considered two different fits having (a) approximately equal pN and ttA branching frac- 
tions, and (b) a larger pN branching fraction. While this choice had little effect on the total width, 
it significantly changed the branching to rjN. 

In order to extract resonance parameters from our global fits, we have generally extrapolated into 
the complex energy plane to search for poles. We have also fitted our energy-dependent and single- 
energy partial-wave amplitudes with Breit-Wigner plus background forms. Here we have chosen to fit 
the partial waves containing the iV(1535) and -/V(1520) resonances in terms of a K-matrix resonance 
form, allowing for 2 poles in the 5*11 partial wave, plus background. This revised parametrization 
was then fitted directly to the data (from 400 to 900 MeV in the pion kinetic energy) in order to 
determine resonance parameters. The remaining partial waves were fixed to values determined from 
a previous global (energy-dependent) fit. This resulted in error estimates more directly tied to the 
data. 

The general form used for the modified (6*11 and -D13) partial waves was 

T = p^l\ (5) 

with giving the phase space for a channel (ttA^, ttA, piV, or r\N\ and with T^. represented in 
terms of a K-matrix as 

= if, (1 + ix,)-' , (6) 

where, in a two-resonance case, we have fitted 

Wi-W W2-W 

The background has been parametrized in terms of the phase space, X^"' = {p^p>y/^K^^ , with k^^ 
elements assumed constant over the limited energy ranges of these fits. The K-matrix pole residues 
were similarly parametrized as =^^^^ with 7' =(/9TY2)^/^. The phase phase factors were 
normalized to unity at the resonance position (Wi). 

Results for the 5*11 and Dis TriV elastic scattering amplitudes are displayed in Fig. |21 Here the 
result of our most recently published fit (FA02) is compared to an updated and improved version 
(G380). The K-matrix fits closely follow the G380 result and have not been plotted. Fig. 0] shows 
the much larger deviations existing between different versions of the i:~p — > ryn amplitudes. 



IV. RESULTS AND DISCUSSION 



A. ?7iV couplings 

Our results from four fits, two with and two without the recent Crystal Ball data, are summarized 
in Tables ITllHIl Listed are the partial widths for the iV(1535) and iV(1520) resonances. While the 
PDG quotes (if a broad and conservative range of about 30 to 50 percent for both the TriV and 
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TjN branching ratios corresponding to the A^(1535), most recent determinations have found the r/N 
fraction to be about 50%, the remaining 50% divided between the ttN and ttttN channels. We have 
similarly found rjN branching fractions exceeding the TriV fraction in all of our fits. The -/V(1535) 
total width, found in the K-matrix fits, differs significantly from our previously published result 53]. 
This is due to the coupled-channel K-matrix form more than a qualitative difference in the partial- 
wave amplitudes. We note that coupled-channel fits have in the past |5^ found the N(1535) width 
to be about half the 200 MeV obtained in single-channel fits to eta photoproduction |57| . 

The extracted A^(1520) rjN branching fraction is very small, as was expected. Penner and Mosel 
found 0.0023±0.0004 for this ratio using an older version of our ttN amplitudes as a representation 
of the ttN elastic scattering database. A somewhat smaller value, 0.0008±0.0001 was found in a 
Mainz analysis of eta photoproduction data Q . We find this quantity to be rather sensitive to model 
details, but our range of values effectively spans the two previous determinations. Fits A and C have 
included the Crystal Ball data, and result in more precisely determined r]N branching fractions, as 
expected. In fits B and D, the Crystal Ball data were excluded. Different pN and ttA branching 
fractions (Fits A and B versus C and D) were obtained through the choice of contributions to the 
background. The background K-matrix contains elements coupling, for example, rjN — s- pN, which 
cannot be measured and are therefore intrinsically mo del- dependent. 

B. 77 A*' scattering length 

As can be seen in Table U previous determinations of the rjN scattering length have produced 
widely varying results. This should not be surprising, as these determinations require the threshold 
behavior of an amplitude that cannot be directly measured. Somewhat surprising to us, however, 
was the relative stability of scattering lengths found from our set of four K-matrix fits, with and 
without the Crystal Ball data, which are shown in Tabic Hvl These results are comparable to those 
found in a fit by Green and Wycech | l48i | , who used a similar K-matrix representation for the S-wave 
amplitudes, and the multi-channel fits of Penner and Mosel 0|- If, however, we determine the r]N 
scattering length directly from our global fit, based on a Chew-Mandlestam K-matrix formalism, a 
very different result is found. This value seems more compatible with the calculation of Ref. d6]. 
As a result, we can confirm previous determinations within similar approaches, but caution that (for 
the real part in particular) the employed model may be more important than improvements in the 
fitted data. The S-wave rjn elastic partial waves, from G380 and Fit A, are displayed in Fig. [S] 

V. CONCLUSIONS AND FUTURE PROSPECTS 

We have explored the model and data dependence of rjN couphngs to the 7V(1535) and iV(1520) 
resonances, and have extracted the ijN scattering length. Our values for these quantities are in 
reasonable agreement with previous determinations. One notable difference in our method has been 
the direct fit to data, rather than to amplitudes. This has allowed a direct comparison of the fits. 

From an experimental point of view, several issues remain to be resolved. The recent Crystal Ball 
measurements of n^p rjn, covering a region from threshold to the peak of the 7V(1535) resonance, 
have suggested a slightly lower mass and width for this state. This could have been verified with 
measurements continuing to higher energies. However, with the Crystal Ball moved to Mainz, this is 
no longer possible. The standard value [l| (547.75±0.12 MeV) of the rj mass has also shifted recently, 
and this naturally effects extrapolations associated with the scattering length determination. The 
linear plot in Fig.|21has taken 547.3 MeV for the eta meson mass. We have allowed the eta mass to 
vary between 547 and 548 MeV, finding very little sensitivity in our fits. We should also note that a 
recent measurement from the GEM collaboration (58j finds a value close to the previous "standard" 
mass of 547.3 McV. 

Given the model dependence found in our determinations, it would be interesting to see the effect 
of the Crystal Ball data in multi-channel fits which include representations of the ttN ttttN data. 
Also of interest would be a re-examination of the A^(1520) i]n branching fraction as extracted from 
eta-photoproduction data. The quality and quantity of data for this reaction has increased since the 
Mainz analysis (3|j . Work on this subject is in progress. 
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TABLE I: r;Af-scattering length overview 



(fm) 


Reference 


AnN (fm) 


Reference 


-0.15 +i 0.22 


[21j 


0.56 +i 0.22 


[21j 


0.20 +i 0.26 


[22] 


0.577 +i 0.216 


[331 


>i 0.24(2) 


[8] 


0.579 +i 0.399 


[24] 


0.25 +i 0.16 


[231 


0.621(40) +i 0.306(34) 


[381 


0.27 +i 0.22 


[5] 


0.68 +i 0.24 


[39] 


0.28 +i 0.19 


[5] 


0.717(30) +i 0.263(25) 


[40] 


0.281 +i 0.360 


[241 


0.734(26) +i 0.269(19) 


[41] 


<0.30 


[25] 


0.75(4) +i 0.27(3) 


[42] 


0.32 +i 0.25 


\M 


>0.75 


[42] 


0.404(117)+i 0.343(58) 


[27] 


0.75 +i 0.27 


[441 


0.42 +i 0.34 


[28] 


0.772(5) +i 0.217(3) 


[45] 


0.42 +i 0.32 


[29] 


0.83 +i 0.35 


[46] 


0.430 +i 0.394 


[24] 


0.87 +i 0.27 


[471 


0.46(9) +i 0.18(3) 


[30] 


0.876(47) +i 0.274(39) 


[27] 


0.476 +i 0.279 


[31J 


0.886(47) +i 0.274(39) 


[27] 


0.476 +i 0.279 


[32] 


0.91(6) +i 0.27(2) 


[48] 


0.487 +i 0.171 


[33] 


0.91(3) +i 0.29(4) 


[49] 


0.51 +i 0.21 


m 


0.968 +i 0.281 


[22 


0.52 +i 0.25 


[35] 


0.980 +i 0.37 


[50] 


0.54 +i 0.49 


[36] 


0.991 +i 0.347 


[2] 


0.55(20) +i 0.30 


[22] 


1.05 +i 0.27 


[42 


0.550 +i 0.300 


[34J 







TABLE II: The present (G380) and previous (FA02 [s^ energy-dependent partial-wave analyses of elastic 
TT^p, charge-exchange (vr'^n), and vr^p ^ ryn (r;n) scattering data, compared to fits A - D from 400 - 
900 MeV. 



Solution 








FA02 


6286/2773 


1920/1100 


635/257 


G380 


5825/2773 


1723/1100 


569/257 


Fit A 


5961/2773 


1684/1100 


539/257 


Fit B 


5935/2773 


1748/1100 


575/257 


Fit C 


6001/2773 


1732/1100 


571/257 


Fit D 


5961/2773 


1839/1100 


582/257 
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TABLE III: Resonance widths (in MeV) and branching fractions. 



Resonance 


Solution 






r,rA 


Fpjv 


rjj/Ftot 


N(1535) 


Fit A 


30±2 


45±3 


15±1 




0.50 




Fit B 


32±3 


45±4 


16±1 




0.48 




Fit C 


39±3 


67±4 


9±2 




0.58 




Fit D 


42±6 


70±10 


11±2 




0.57 


N(1520) 


Fit A 


68±1 


0.12±0.03 


19±5 


19±5 


0.0012 




Fit B 


68±1 


0.17±0.12 


19±6 


19±6 


0.0016 




Fit C 


67±1 


0.08±0.03 


14±4 


24±4 


0.0008 




Fit D 


(.)7±i 


0.09±0.07 


il±5 


2i±5 


O.OOQ!) 



TABLE IV: scattering lengths from K-matrix fits (resonance plus background, see text) and the global 
energy-dependent fit (G380). 



Solution 


Scattering Length (fm) 


Fit A 


1.14 + i 0.31 


Fit B 


1.10 + i 0.30 


Fit C 


1.12 + i 0.39 


Fit D 


1.03 + i 0.41 


G380 


0.41 + i 0.56 
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100 





80 


I/} 




60 


C 


40 


XI 




20 









200 




160 


m 


■s 
-i- 


120 


C 


80 


Xj 


o 
T) 


40 









300 


Sh 


250 


t/] 






o r\ o 

2UU 


3 






150 


c 




XI 


100 


b 




X) 


50 









300 




250 


m 






200 


3 




150 


G 




X) 


100 


b 


50 


X3 



561 MeV 



(a) 



i i 



100 



-0.5 



0.0 

COS0 



0.5 



1.0 




578 MeV 



(c) 



1.0 



-0.5 



0.0 

COS0 



0.5 



1.0 




600 MeV 



(e) 



-0.5 



0.0 

COS0 



0.5 
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FIG. 1: (a)-(g) Differential cross sections for tt p — > ryn at seven incident vr energies. The uncertainties are 
statistical only. For the total cross sections (h), we have combined statistical and systematic uncertainties 
in quadrature. FA02 (E913/E914 data not included), G380, and Fit A shown as solid, dash-dotted, 
and dotted lines , re spectively. Experimental data are from ;7*| (filled circles), (open circles), [TH (open 
triangles), and (open squares) measurements. Other previous measurements (for references see SAID 
database 1541) shown as asterisks. 
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FIG. 3: (a) Sn and (b) D13 partial amplitudes for ttN elastic scattering. Solid (dashed) curves give the real 
(imaginary) parts of amplitudes corresponding to the predictions of solution FA02 [ssll (E913/E914 data not 
included). Single-energy solutions associated with FA02 are plotted as filled and open circles. Dash-dotted 
(dotted) curves show the real (imaginary) parts of amplitudes corresponding to G380. Differences between 
G380 and Fit A are not significant. All amplitudes are dimensionless. Vertical arrows indicate Wr and 
horizontal bars show full T/2 and partial widths for F-^jv associated with the FA02 results. 




FIG. 4: (a) Su and (b) Dra partial amplitude for n^p rjn. Dash-dotted (dotted) curves show the real 
(imaginary) parts of amplitudes corresponding to G380. Solid (short-dash-dotted) lines represent the real 
(imaginary) parts of amplitudes corresponding to the Fit A. All amplitudes are dimensionless. 
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FIG. 5: -p^ dependence of the Sn amplitude for the reaction rjn rjn. Dash-dotted (dotted) curves give 
the real (imaginary) parts of amplitudes corresponding to the solution G380. Solid (short-dash-dotted) lines 
represent the real (imaginary) parts of amplitudes Fit A. All amplitudes are dimensionless. 



